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Abstract. We present a method to measure the critical temperature of the island of a 
superconducting single electron transistor. The method is based on a sharp change in the slope of 
the zero-bias conductance as a function of temperature. We have used this method to determine 
the superconducting phase transition temperature of the Nb island of an superconducting single 
electron transistor with Al leads. We obtain T^ b as high as 8.5 K and gap energies up to 
ANb — 1-45 meV. By looking at the zero bias conductance as a function of magnetic field instead 
of temperature, also the critical field of the island can be determined. Using the orthodox theory, 
we have performed extensive numerical simulations of charge transport properties in the SET 
at temperatures comparable to the gap, which match very well the data, therefore providing a 
solid theoretical basis for our method. 
PACS: 73.23.Hk,73.40.Gk,74.50.+r 



The single-electron transistor is a device with remarkable properties, which have been 
intensively studied over the last almost two decades [1]. The first application of such structures 
have been in the area of electrometry [2], but recently new developments have promoted the 
superconducting version of this transistor as one of the main contenders for realizing the so-called 
charge quantum bits [3]. These devices could become the building blocks of future quantum 
computers, in architectures that allow exchange of quantum information between them via 
transmission lines [4]. 

These new areas of interest have created a demand for high performance nanofabrication 
techniques and there has been a strong motivation to develop a technique for a reliable fabrication 
of a small niobium-based Josephson junction. Niobium (Nb) would provide more reliable 
performance in many of these devices due to its large superconducting gap ANb ~ 1-5 meV 
(in bulk) as compared to aluminium, Aai ~ 0.2 meV, which has been the material of choice in 
nanofabrication for many years due to its easier processability. A larger superconducting gap 
would provide a better suppression of the undesirable quasiparticle tunnelling and also would 
yield larger Josephson coupling energy Ej, in ultrasmall junctions. 

The conventional e-beam shadow evaporation (self- alignment) technique, has been applied 
successfully for soft metals like Al, Cu and Pb. Yet, it is known to be difficult to apply this 
method for refractory metals like Nb due to the poor quality of the produced films, resulting in a 
critical temperature T c well below the bulk value. To circumvent this problem, several fabrication 
methods have been proposed, and different groups have been able to obtain good quality small- 



capacity single Nb-Al and Nb-Nb junctions [5] and Nb-based single electron transistors [6]. 

In all these devices, since there is no easy direct electrical access to the island, the typical 
way to determine the critical temperature is to produce separate samples with films of the same 
thickness, and measure the critical temperature of those. However, the critical temperature of 
the Nb island can be quite different from that of a film, since the island's dimensions are small, 
comparable with the superconducting coherence length (^Nb = 38 nm in the bulk). Of course, 
determination of the gap at low temperatures from tunneling data provides as well an indirect 
measure of the critical temperature via the BCS relation: but again, for lower-dimensional 
superconducting objects the corrections to this law could be significant. 

In this article we present a method to determine directly the critical temperature, and also the 
critical magnetic field, of a Nb island in superconducting Al/AlOx/Nb/AlOx/Al-SETs fabricated 
with the standard self-alignment lithography process using a slightly modified recipe. 

The fabrication technology for our samples has been described elsewhere [5, 7]. The 
measurements were carried out in a small dilution refrigerator (Nanoway, PDR50) with the 
base temperature below 100 mK and well-filtered electrical lines. 

In Fig. 1(a) we show an SEM picture of one of the samples and its IV characteristic. In the 
SEM image, brighter lines are niobium and darker lines are aluminium. Besides the existence 
of a quasiparticle voltage threshold and the Josepshon current, the IV displays a number of 
features, such as Cooper-pair resonances and low-bias excess currents which have been analyzed 
elsewhere [7]. 
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Figure 1. (a) SEM image of one of the measured samples with inset showing the schematics 
used for measuring the SET. (b) The I-V characteristics and dl / dV-cuive measured from a 
sample with Eq ks 35 fxeV and Ej ks 92 fxeV. 

The charging energies of the samples were derived from the normal state conductance curve 
measured at 4.2 K with magnetic field of B ~ 5 T. The obtained charging energies varied 
between Eq = e 2 /2C sa 34 — 78 yueV. The Josephson coupling energy was found from the 
Ambegaogar-Baratoff formula for two different superconductors, which yielded Ej ~ 41 — 123 
^eV. 

In Fig. 1(b) the I-V characteristics and a dl /dV -curve measured at the temperature T ~ 200 
mK are shown from one of the samples with Eq ~ 35 //eV and Ej ~ 92 fieV. The gap - with a 
width of 4(Aai+ ANb) /e - is clearly visible and the maxima in the dl /dV-curve yield Ap>jb ~ 1-45 
meV. Here we have assumed Aai ~ 0.2 meV. At bias voltages below V ~ 0.9 mV we also see 
(Fig. 1(b), the IV curve) a series of gate-dependent, equally-spaced peaks, which appear only 
when the aluminium leads become superconducting. They correspond to Cooper pair resonances 
occuring when charges are transported in the whole circuit [7]. 

Experimentally, our method consists of determining the critical temperature T^ h as well 
the upper critical field # c N 2 b fr om the zero bias conductance of the sample as a function of 




temperature or magnetic field, respectively. An example is shown in Fig. 2 where the sudden 
transition from the strong temperature dependence determined by the reduction in A^iT), to 
a much weaker temperature dependence (due only to Coulomb blockade) is clearly visible at 
the phase transition. The critical temperatures obtained were 7.8, 8.1 and 8.5 K in the three 
best samples and above 7.5 K in all others. This is already very close to T c Nb w 9 K in bulk 
Nb. The obtained critical fields were between # c N 2 b « 2.5 - 4.5 T. In this case the transition 
to the normal state happens gradually between the thermal critical field and the upper 
critical field where superconductivity is completely suppressed (thin, disordered Nb films 
are type-II superconductors), as seen on the left inset of Fig. 2(a). 
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Figure 2. (a) Example of determining the critical temperature T c in one of the samples. Circles 
are the measured zero bias conductances obtained from the experimental dl /dV -curves shown 
in the lower-right inset. The dashed line is a linear fit to the data below T c . The slope is 
determined by the reduction of ANb(^) as the temperature is increased. The dotted line is a 
linear fit to the data above T c and is rather flat due to the much weaker temperature dependence 
of the Coulomb blockade. In the left inset we present similar data for magnetic field instead of 
temperature, a graph which we used for determining the upper critical field of Nb island, 
(b) Data points (circles) and numerical calculations for charging energy, Eq = 35 //eV. The 
continuous line corresponds to junctions without excess current, that is a = V = 0; the dotted 
line was calculated with a = 0.38 meV, T = 37.5 /ueV, and a best fit around the critical 
temperature was obtained (dashed line) for a = 0.38 meV, T = 161.25 fieV. 

To model or data, we have used the standard orthodox theory [1] for charge transport in NISIN 
single-electron transistors (at the temperatures of interest, around the critical temperature of 
Nb, the Al leads are normal). We calculate numerically the tunneling probabilities in each 
junction for a given superconducting gap, and to find the currents and conductances at a given 
bias voltage we solve the corresponding master equation for charge transport across the structure 
for each given temperature. From the low-bias Cooper pair resonance data (which depend on 
the charging energy) we extract, as in [7], Eq = 35 /xeV. For the gap, which is used to calculate 
the tunneling probabilities, we use the standard BCS gap equation at finite temperature 

1 7 tanh(/V^ + A(T)V2) 

N{0)V ~ J € y/S + A(T)* 

where N(0) is the density of states at the Fermi level, V is the BCS interaction, (3 = 1/ksT and 



the Debye energy is taken much larger than the gap. Since A(T = Tq) = 0, we can write 
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an integral equation which is used to determine numerically the gap value at a given temperature 
T. 

Also, when calculating the tunneling probabilities, in order to take into account the excess 
subgap currents present in Nb-based junctions, we have used a life-time broadening V of the 
quasiparticle energies, resulting in a density of states [7] 
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J o 



Re 



E-iT 



(1) 



where V is the probability density associated with a certain value of the gap, where the function 
V is a Gaussian of standard distribution a. The value of a was determined from this sample 
using the low-bias IV plots, as described in [7]: its value is set by the fabrication process and 
does not depend on temperature. The life-time broadening T depends in general on temperature: 
its value at T = 95.7 mK was detemined in [7] as T = 37.5 peV. We assume that this dependence 
is mild near the phase transition so that T can be considered constant: we find, as a best fit to 
the data, T = 161.25 fieV. We note also that at 5-6 K this value does not give a good fit, and 
a better match for the data would require indeed a lower value of T, corresponding to a lower 
temperature. Fig. 2(b) summarizes these findings. 

We described a new method for determining the critical temperature of a Nb island which 
forms the active element in a superconducting single electron transistor. The method is based 
on the experimental observation that, as the system is cycled through the phase transition, the 
slope of the zero-bias conductance as a function of temperature changes abruptly from a finite 
value (in the superconducting state) to almost zero (in the normal state). 
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